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Dirac electrons in solids show characteristic physical properties due to their linear dispersion relation and two-band 
nature. Although the transport phenomena of Dirac electrons in a normal state have intensively been studied, the transport 
phenomena in a superconducting state have not been fully understood. In particular, it is not clear whether Dirac electrons 
in a superconducting state show Meissner effect (ME), since a diamagnetic term of a current operator is absent as a result 
of the linear dispersion. We investigate the ME of three dimensional massive Dirac electrons in a superconducting state 
on the basis of Kubo formula, and clarify that Meissner kernel becomes finite by use of the inter-band contribution. 
This mechanism of the ME for Dirac electrons is completely different from that for the electrons in usual metals. Our 
result shows that the Meissner kernel remains finite even when the superconducting gap vanishes. This is an unavoidable 
problem in the Dirac electron system as reported in the previous works. Thus, we use a prescription in which we subtract 
the normal state contribution. In order to justify this prescription, we develop a specific model where the Meissner kernel 
is obtained by the prescription. We also derive the result for the electron gas by taking the non-relativistic limit of Dirac 
Hamiltonian, and clarify that the diamagnetic term of the Meissner kernel can be regarded as the inter-band contribution 
between electrons and positrons in terms of the Dirac model. 


1. Introduction 

Dirac electron systems have been attracting a great inter¬ 
est in recent years. One of the characteristic features of the 
Dirac electron systems is that its motion is described by Dirac 
equation instead of Schrodinger equation. Although the orig¬ 
inal Dirac equation describes the motion of relativistic (i.e., 
high-energy) electrons, 1 * it is known that the electronic states 
of several materials with linear dispersion relation near the 
Fermi energy can be described by a low-energy effective 
model which has the same structure as Dirac equation. Exam¬ 
ples of such materials are graphene, 2 * a-(BEDT-TTF) 2 l 3 , 3 * 
Bi, 4 * Ca 3 PbO, 5 * iron pnictide, 6,7 * Na 3 Bi, 8 * Cd 3 As 2 , 9 * and so 
on. It has been reported that the physical properties of these 
materials are qualitatively different from those in the usual 
metals which are described by Schrodinger equation and 
Fermi liquid theory. Hence, it is important to investigate the 
fundamental properties of the Dirac electrons. 

Up to now, many works have been carried out to inves¬ 
tigate the normal-state properties of Dirac electron systems, 
such as Hall conductivity, 10 * magnetoresistance, 11 ’ 12 * orbital 
magnetism, 13-15 * Nernst coefficeint, 16 * and spin Hall conduc¬ 
tivity. 17 - 181 In these studies, it has been pointed out that the 
inter-band effect plays important roles, since the conduction 
and valence bands are close to each other. For instance, the 
giant orbital magnetism in Bi is due to the inter-band effect, 
and can not be explained by Fandau-Periels’s theory 19 - 20 * of 
diamagnetism in which the inter-band is neglected. 

On the other hand, the transport properties of Dirac elec¬ 
trons in a superconducting state have not been understood 
completely. In particular, it is not trivial whether the Dirac 
electrons in the superconducting state show the Meissner ef¬ 
fect (ME), since the current operator does not have the dia¬ 
magnetic term; 21 * for usual metals with a parabolic dispersion 
relation (i.e., in electron gas), the ME arises from the diamag¬ 
netic term of the current operator. Note that the diamagnetic 


current term appears when the Hamiltonian has a kinetic en¬ 
ergy proportional to k 2 with k being the momentum. 22 * It is 
therefore necessary to discuss the ME of Dirac electron in a 
superconducting state on the basis of the treatment which cor¬ 
rectly includes the inter-band effect. 

Although there have been many works on the supercon¬ 
ductivity in the Dirac electron systems, there are only a few 
studies on the ME. 23,241 In Ref. 23, the coexistence of the 
nodal charge density wave (CDW) and the superconductiv¬ 
ity is studied. Since the excitation spectrum of the mean-field 
CDW state has Dirac-like dispersion near the Fermi level, the 
situation can be regarded as a superconductivity of a Dirac 
electron system. However, their result shows that the Meiss¬ 
ner kernel at zero temperature vanishes when the Fermi level 
is on the Dirac point although the superconducting order pa¬ 
rameter A is nonzero. This is rather unphysical. As we show 
in the present paper, the inter-band contribution gives a finite 
Meissner kernel even when the density of states at the Fermi 
level is zero. This indicates that the inter-band contribution is 
not taken into account correctly in Ref. 23. 

In Ref. 24, the Meissner kernel of the superconducting state 
of graphene is calculated. However, the obtained result shows 
that the ME remains even in the normal state. This is appar¬ 
ently unphysical. The authors of Ref. 24 suggests that this 
problem can be resolved by subtracting the normal state con¬ 
tribution of the Meissner kernel. However, they do not obtain 
the Meissner kernel which remains after subtraction. 

In this paper, we discuss the ME of a 4 x 4 massive Dirac 
electron system in the three-dimensional space assuming a s- 
wave Cooper pairing. The derivation of the Meissner kernel 
is based on Kubo formula. We reveal that the inter-band ef¬ 
fect plays an important role and that the ME appears as an 
inter-band contribution in spite of the absence of the diamag¬ 
netic term of the current operator. Our results also show that 
the ME remains in the normal state as in Ref. 24. This prob¬ 
lem will be unavoidable as far as we use the unbounded Dirac 
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dispersion. Therefore, we use the prescription to subtract the 
normal state contribution of the Meissner kernel as in Ref. 
24. We discuss the obtained Meissner kernel as a function of 
chemical potential. In order to justify this prescription, we de¬ 
velop in Appendix a model which is an extension of the Dirac 
Hamiltonian. In this specific model, we show that the Meiss¬ 
ner kernel vanishes in the normal state, and that the correct 
Meissner kernel in the superconducting state is obtained by 
the prescription. Although we use a specific model, we expect 
that this prescription is reasonable. 

As discussed above, the inter-band effects in Dirac electron 
systems have been studied in the normal state. On the other 
hand, it has not been recognized so much that the inter-band 
effects play important roles in the response functions in the 
superconducting state. In this paper, we show that the inter¬ 
band effect is essential to the Meissner kernel, on the basis of 
the Kubo formula. 

This paper is organized as follows. In the next section, we 
introduce a 4 x 4 massive Dirac Hamiltonian in the three- 
dimensional space, and the definition of the current operator 
in that model. We also show the treatment of the supercon¬ 
ducting order parameter in the mean field approximation. In 
Sect. 3, we show the explicit form of the Meissner kernel in 
the present model by use of Kubo formula, and give the re¬ 
sult of the numerical analysis. Then we discuss how Dirac 
electrons become the Meissner state. It will be shown that the 
inter-band contribution plays an important role in obtaining 
the finite Meissner kernel. We also mention the relation be¬ 
tween the Dirac electron and the non-relativistic electron gas 
by considering the large band gap limit case. We will show 
that our theory can reproduce the well-known results in the 
non-relativistic electron gas, and that the origin of what we 
call the paramagnetic and diamagnetic terms originates from 
the intra and inter band term, respectively. Finally, the brief 
summary is given in Sect. 4. In Appendix, we develop a spe¬ 
cific model in which the normal state Meissner kernel van¬ 
ishes. Part of the present work has been published before. 251 


2. Formulation 

2.1 Hamiltonian 

We consider the following 4x4 massive Dirac Hamiltonian 
in the three-dimensional space: 10 ' 


H 0 {k) 



MI 

—ivk ■ cr 


ivk ■ cr \ „ 
—MI ) ° kl 


(1) 


where I is the 2x2 unit matrix and cr = (<j T , o y . <7 : ) are Pauli 
matrices. M is the band gap at the expanding center in the 
Brillouin zone, and v is the Fermi velocity. The basis used in 
(1) is 4 = Where 1,2 denote the 

orbital, and f, j denote labels of the time reversal pair. 

Diagonalizing this Hamiltonian, we obtain the following 
Hamiltonian, 

»»(*>=«:( T <2) 



Fig. 1. (Color online) The dispersion relation of the 4x4 Dirac Hamilto¬ 
nian in normal state. Red line corresponds to the upper band with index +, 
and blue line corresponds to the lower band with index —. Each band has 
two-fold degeneracy labelled by -ft and JJ-. 


1]. The unitary matrix which is defined as 4 
given by 


U(k) = 


X{k)I —iY{k) cr 
—iY(k) ■ (T X(k)I 


= U ( k)c k is 


( 3 ) 


whereX(t)= /and Y(k)= 2e{k)(e \ k)+M) vk. 


2.2 Current operator 

In this subsection, we discuss the current operator in the ab¬ 
sence and presence of the vector potential. Without the vector 
potential. A, the current operator in the momentum space is 
given by 


j{q)=-eY, d \- q d kHo(k)c k 

k 

=- e Yc : ( 0 iwr 

Z. k-q ^ _i vo - 0 


4 


(4) 

k 


where e > 0 is the absolute value of the charge of an electron. 
The vector potential A is introduced in the Hamiltonian by 
replacing k by k + eA. However, we can readily see that the 
current operator in Eq. (4) does not change even in the pres¬ 
ence of A, since Ho(k) contains only linear terms with respect 
to k. 

This property of the current operator is essentially differ¬ 
ent from that in the electron gas in which the current operator 
is given by j(q) ='Lk.c c l-q, a 7 ,i k + eA ) c k,a where the last 
term is called as the diamagnetic term. Apparently, the cur¬ 
rent operator of the Dirac electron systems does not have dia¬ 
magnetic current. As we will show shortly, the absence of the 
diamagnetic term is crucial for discussing the mechanism of 
the ME in Dirac electron in superconducting state. 


where e(k) = \/M 2 + (yk) 2 , and a k represents the band-basis 
which is expressed as a k = {a ki+ ^,a k&A ,a k -^,a k ,^) T . 
Here, the indices + and — denote the upper and the lower 
bands respectively, and jj and JJ. denote psudo-spins corre¬ 
sponding to the two-fold degeneracy of each band [See Fig. 


2.3 Superconductivity 

In this subsection, we introduce superconductivity in the 
Dirac electron system. We assume the following attractive in- 
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teraction, 

Hml = —V 52 12 a l,ri,f a -k,i]4 a -k',r]' ( 5 ) 

|*|,|* / |<* c ?7,17 , =± 

where k c is the cutoff momentum corresponding to the range 
of energy in which the attractive interaction works. When 
vk c 3> M, a superconducting state is realized even if the chem¬ 
ical potential is located in the band gap (i.e., there is no den¬ 
sity of states at the Fermi level) for sufficiently large V, since 
electrons with k<k c can still contribute to the Cooper pairing. 

We define the superconducting order parameter in the s- 
wave symmetry (see Fig. 2(a)), 

A± = V 52 ( a -k',±,ij, a k',±, ft)- (6) 

\k'\<k c 

By applying the mean field approximation to the Hamiltonian, 
H = Hq + //i nt , we obtain the BCS Hamiltonian, 

tfBCS = 52 E 12 ^±( k ) a lr 1 Z a k.riX 

k J?=±£=M 

52 12 (7) 

k i)=± 





Fig. 2. (Color online) (a) Schematic picture of the cooper pairing, and (b) 
energy eigenvalues for Bogoliubov quasi-particles. 


each other. 


where A (k) = (A+ + A_)0(k c — |£|) = A &(k c — |&|) with 
0 (jc) being the step function, and t^±(k) = ±e(k) — fi. 

Diagonalizing //iic:;s - we obtain the the excitation energy 
E±(k) = y/^±(k) 2 + A 2 (k) [see Fig. 2(b)], Then, thermal 
Green’s function, 

[y{k,ico„)Ua 2 = -J q P dre i(O " x (T T [a Kai (T)al a2 (0)}), (8) 
and the anomalous Green’s functions, 


rP 

[^{k,i(O n )]a l a 2 = — J dx e lCOl,l: (T z [ah ai (T)a_fc iO [ 2 (0)]), 

(9a) 

(k,iCO n )\ ai a 2 = -J^ dxe l ^(T T [al kai (x)a^ a2 {0)]), 

(9b) 

are obtained in the form of 4 x 4 matrices as follows. (Note 
that O'! and a x are the sets of band and pseudo spin indices.) 


&(k,i(On) = 


y+{k,i(D n )i 

0 


Sf-( kja> n )l J ’ 


and 

cP(k,i(D n ) = — ^(k,iC0n) 

ic y ,^ + (k.i(o n ) 

0 


where 


and 


y±{k,i(D n ) = - 


ui(k) 


io y ^-(k,i(On) 


v±(k) 


&±(k,i(On ) = -u±(k)v±{k)(- 


ico n —E±(k) ico n + E±(k )’ 

1 1 


00 ) 


( 11 ) 


( 12 ) 


i(D n — E±(k) iatn + E±’(k) 

’ 03 ) 


with ul(k) = \[l + v|(*) = i[l - jf^], and 

u±(k)v±(k ) = 2 E±(k) • ^ should be noted that the upper and 
lower parts of these Green’s functions are decoupled with 


3. Meissner effect in Dirac electron systems 

3.1 Small band gap case: Dirac electrons in solids 

Using the formulations derived in the previous section, 
we calculate Meissner kernel, K xx , through Kubo formula. 

Since there is no diamagnetic term of the current operator, 

(S) 

K xx {(\.ioti) is given only by the current-current correlation 
function, 

K^\ qi m x ) = - f P dTe im ^(T T (j x (q,T)j x (-q,0))), (14) 
Jo 

where co x = (A = 0,±1,±2,---) denotes a bosonic 
Matubara frequency and T an imaginary time. The Meissner 
kernel for real frequency is given by the analytic continuation 

^(q,®) = K^\q,icox)\ ia>k ^w 

Applying the Bloch-De Dominies theorem to Eq. (14), we 
obtain the following expression of K xx (q,ico x ): 

K xx \q,iC 0 X ) 

=T 52 Tr [&(k- q,<On - (0x)J x {k,q)&{k,(0n)}x{k- q,-q) 
k,m„ 

-^\k-q,(0 n -(0 X )J x {k,q)^ r (k,(0n)j^{-k,-q)], (15) 

where J x (k„q) is given by J x (k,q) = [U(k - q)j x U ^(*)], 
and co n = ( 2 ”+ 1 ) ?r ^ w = 0,±1,±2,-") denotes a fermionic 
Matubara frequency. It should be noted that the limit co —> 0 
has to be taken before the limit q —tO. 

By taking the trace in Eq. (15), we obtain the Meissner ker¬ 
nel which consists of two parts as 

K%\q, 0) = KT XS \q, 0) +K^ S \q, 0). (16) 

^.inua.fS} (q o) comes f rom ^e term in which the two Green’s 
functions in Eq. (15), (S?Sf or .'¥ ■'¥) have the same band in- 
deces (+,+) or and K™ ter '^ S \q,0) from the Green’s 

functions with the opposite band indeces, i.e., (+,—) or 
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Fig. 3. (Color online) M dependence of a// (0,0), (0,0), and K- a ( 0,0) for (a) A — 0.1M, (b) A = 0.1M, and (c) A (k) = 0.7M. M dependence of K xx is 

in the inset. 


The explicit forms of these are obtained as 

^r (s) ( ,,o 


*£"(*>«) = 


+ 


e(k)e(k - q) 

e(k)e(k — q) — M 2 
e(k)e(k — q) 




v 2 [k x (k x q x ) ky(ky (/>) k z (k z q z )\ 


e(k)e{k - (/) 
e(k)e(k — q) +M 2 


e(k)e(k - q) 

In Eqs. (17) and (18), B ’s are coherence factors, 


B 


ter(i) (jti q) = i + Z + (k-q)t;-(k)-A 2 (k) , 


E+(k-q)E-(k ) 


(Si 

When we take the limit of q — ► 0, we obtain A+7(0,0) = 
Ar”“"' S) (0,0)+ C ter ’ (s) (0,0) with 


k rj=± 

.. rpintra(i) , . , /(^l] (*)) “ “ «)) 

„intra(ii),, , f{^i\ W) + /(^T) (A — ?)) — 1, 

11 £„(*-*)+£„(*) (17) 

and 

4r |s, («,o)=£f ££ s r(M) 

^ ^ t?—i 

X (t,,)/(E-,(t))-/(gM*-<» 

1 V E z] {k~q) —E- rl (k) 

-B inl ^(k,q) nE ~l {k ^ +/ S* 7/' )} ~ 1 ]■ (18) 

E n {k- q)+E- n (k) 

Here, /(x) is the Fermi distribution fucntion, and g l ^ ra {k,q), 
g^ el (k,q) are given by 

f [k x (k x qx ) — ky(ky ) k z (k z <+)] 


C tra ’ (s) (0,0) = 2 e V 




and 


/r 1 

A, 


inter, (S) 


r W + (*)) , df{E-{k)) 

dE-(k) J 

( 21 ) 


v 2 k 2 

e 2 (k) 


(0,0) = -2eV£ 1- 


x[ 1 + 


(19a) 


(19b) 


S + m_(k)-A\k) \ f(E-(k))-f(E + (k)) 
E+{k)E-{k) J E+(k)-E-(k) 

_ %+{k)^{k)-A 2 {k) \ f{E-(k))+f(E+(k))- 1 
E+{k)E-{k) J £+(*)+£_(*) J ' 

( 22 ) 

It should be noted that for T = 0, K™ tia ’ (0,0) is 0 since 
the spectrum has a finite gap at the chemical potential. Hence, 
Kxx only consists of the inter-band contribution. This is in 
sharp contrast to the conventional s-wave BCS case, in which 
the Meissner kernel from the paramagnetic current vanishes at 
T = 0, but the Meissner kernel from the diamagnetic current 
remains which leads to the ME. In contrast, in the Dirac case, 
the diamagnetic current does not exist from the beginning, but 
instead the inter-band contribution exists, which leads to the 
ME. 

However, there is an interesting problem appearing from 
the inter-band contribution. Even if we put A (k) = 0, i.e., in 
the normal state, the kernel (0,0) does not vanish. The 
explicit form is 

v 2 ki r df(^ + (k)) , dftt-m 


fd" tra,(N) (o, o) =2 e V£ 


b; „< 0 (m) = 1+ MA_«|W. (20 a, 


T £2 ( k ) 


d^k) 


d^(k) J 


e 2 v 2 (m 2 -M 2 ) 3 / 2 ©(|m|-M) 


= on,, 


3 7l 2 


and 


^ ter ’ (N) (o,o) 


(20c) 


=-2*ve 1- 


v 2 k 2 

e 2 (k) 


|M|v 3 


/(£_(*))-/(&.(*)) 


(23) 


e(k) 


ninter(ii) / i. \ i £+(* ~ q)$- (*) ~ A 2 {k) 

B {k ' q) ~ l - E + {k — q)E_(k) ‘ (2 ° d) 


e 2 v 2 
37r 2 


A 3 (m 2 -M 2 ) 3 / 2 0(|m|-M) 


e(A) 


|M|v 3 


(24) 
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hence 

k£P (0,0) = £;£ tra ’ (N) (0,0)+ tfS ter,(N) (o, 0) 


eV r A 3 - 
3n 2 e(A) 


(25) 


Here we introduce the ultraviolet cut-off momentum A. This 
finite Meissner kernel in the normal state is due to the linear 
dispersion and the unboundedness of the spectra in the Dirac 
electron systems. A similar problem was pointed out before in 
graphene. 24) In the usual electron gas model, K'^' 1 = 0 holds 
since the paramagnetic term exactly cancels with the diamag¬ 
netic term. However, in the case of Dirac electron systems, 
there is no “counter contribution” by which the normal state 
Meissner kernel vanishes. Therefore, in the following, we use 
the prescription to subtract the K^' 1 as suggested in Ref. 24, 
i.e., we calculate the Meissner kernel in the superconducting 
state by 

K a = K^-K^\ (26) 


It should be noted that only the contribution from |£| < k c , 
where A (k) is finite, gives the finite value in Eq. (26), since 
the contribution from k > k c in the normal and supercon¬ 
ducting states are exactly same and cancel with each other. 
In order to justify the prescription Eq. (26), we develop in 
Appendix a model which is an extension of the Dirac Hamil¬ 
tonian and show that the Meissner kernel in the normal state 
vanishes. In this specific model, we also show that the cor¬ 
rect Meissner kernel in the superconducting state is obtained 
by the prescription of (26). We should note that it is not clear 
whether this prescription can be applied to general Dirac elec¬ 
tron systems, although it is justified in the specific model in 
Appendix. However, it is physically required that the Meiss¬ 
ner kernel should vanishes in the normal state and we expect 
that the prescription (26) is reasonable. 

By using Eq. (26), we numerically calculate Xxt(0,0) at 
T = 0. In the numerical calculation, we set k c = lOM/v. 

Figures 3 shows p dependence of ,K& \ and K xx for 
a fixed A, (a) A = 0.1 M, (b) A = 0.5 M , and (c) A = 0.7M 
respectively. (The insets show K xx in the extended scale.) 
We can see that although and 1 are both negative, 
K xx ( 0,0) is always positive and finite, thus we obtain the 
Meissner state. 

We can also see in Figs. 3, that K „ strongly depends on A, 
in ~M < p < M. This is another sharp contrast to the con¬ 
ventional BCS case in which K xx does not depend on A and 
depends only on the carrier density. 


3.2 Large band gap case: non-relativistic electron 

Up to now, we consider the case in which the band gap 
is very small (~ 10 meV). In this subsection, we consider 
the opposite case, i.e., the large band gap case. In particular, 
it is important to test whether the present theory reproduces 
the results for non-relativistic theory. To access to the non- 
relativistic limit, we replace M to me 2 , and v to c, with m being 
the mass of an electron and c the speed of light. Then, we can 
regard the upper band as the real electron, and the lower band 
as the positron (see Fig. 4). We further assume that the chem¬ 
ical potential is located in the upper band, and that p — me 2 is 
much smaller than me 2 . Then we obtain the Meissner kernel 


in the normal state as 


^intra,(N) 

&xx 


( 0 , 0 ) 


and 


^inter, (N) 
&xx 


( 0 , 0 ) 


e 2 c 2 (p 2 — m 2 c 4 ) 3 / 2 e 2 n 
3 K 2 pc 2 m 

e 2 c 2 [ A 3 {p 2 -mV) 3 / 2 
37T 2 e(A) pc 3 


(27) 


e 2 n e 2 c 2 A 3 
m 37r 2 e(A)’ 


(28) 


where kp = \J2m(p — me 2 ) is the Fermi momentum, and n = 
is the electron density. Here we use the approximantion 

p — me 2 <C me 2 . 

Comparing Eqs. (27) and (28) with the results for the elec¬ 
tron gas, 22 - 1 we can see that K™ tm ^ is exactly equal to the 
paramagnetic term in the electron gas and that k'( l ( lel ' ! ;V :i is ex¬ 
actly equal to the diamagnetic term apart from the divergent 
part (—Note that this divergent term comes from 
the region far away from the Fermi surface, and that this is 
exactly what we deal with by the subtraction in order to avoid 
the unphysical result in the last subsection. We think that this 
kind of divergence is inevitable within the Dirac theory. Eqs. 
(27) and (28) indicate, first, that we can reproduce the result 
of the electron gas starting from the Dirac model; second, that 
the diamagnetic term of the kernel in the electron gas origi¬ 
nates from the inter-band contribution between positrons and 
electrons. 

Now we turn to the superconducting state. We assume that 
the range of the energy in which the attractive interaction 
works is much smaller than me 2 or p — me 2 as shown in Fig. 4. 
Thus, the pairing occurs in the very vicinity of the Fermi level 
on the upper band [see the yellow shade in Fig. 4]. In this case, 
kT' {S) becomes 0 as we have discussed in the last subsec¬ 
tion. On the other hand, is almost same as ^-“ ter >( N ), 

since the states below the Fermi level do not change by the 
pairing except for the vicinity of the Fermi surface. There¬ 
fore, we obtain the Meissner kernel in the superconducting 
state as 


K 


(S) 


jointer, (N) 
&xx 


2 

e n 
m 


(29) 


which is completely consistent with the result of the super- 



Fig. 4. (Color online) Dispersion relation in the non-relativistic limit. Yel¬ 
low shade denotes the range where the attractive interaction works. 
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conducting state in electron gas. 

4. Summary 

In summary, we have investigated the ME of three- 
dimensional massive Dirac electron in superconducting state 
on the basis of Kubo formula. Although the diamagnetic term 
of the current operator is absent in Dirac electron system, the 
Meissner kernel finite for any value of the chemical potential, 
since the inter-band contribution remains finite. This inter¬ 
band mechanism of the ME of Dirac electron is essentially 
different from that of usual metals. 

In the Dirac electron system, there is an unavoidable prob¬ 
lem that the Meissner kernel remains finite even in the normal 
state. In order to obtain the Meissner kernel in the supercon¬ 
ducting state, we use a prescription of Eq. (26). Although we 
have not proved that this prescription is completely justified, 
we develop a specific model in which the Meissner kernel 
vanishes in the normal state and the prescription is reason¬ 
able. 

We also have studied the large band-gap limit to discuss the 
non-relativistic case. We clarify the relation that the paramag¬ 
netic term of the kernel originates from the intra-band term 
and the diamagnetic term from the inter-band term. 

TM is grateful to Y. Masaki for fruitful discussions and 
comments. This work was supported by Grants-in-Aid for 
Scientific Research (A) (No.24244053). 


Appendix: Meissner effect in Dirac electron systems 
with a quadratic term 


In this appendix, we perform a simple model calculation 
in which the Meissner kernel in the normal state vanishes. In 
order to avoid the unboundedness of the Dirac dispersion, we 
add a kinetic energy term |k 2 in the Hamiltonian instead of 
introducing the Brillouin zone. 


H 0 {k) 



(M+fk 2 )/ 
—i vk ■ cr 


\vk ■ cr 

(-M+fk 2 )/ 



(AT) 


where parameter a is assumed to be very small. The energy 
dispersion becomes e(k) + fk 2 , and a new Fermi surface ap¬ 
pears as shown in Fig. A-l. The matrix form of the current 
operator, Eq. (4) changes as 


j { q ) = - e E ( 


k q 


a(k + eA)I 
—ivcr 


ivcr 

a(k + eA)I 


Ck- 

(A-2) 


Note that there appears a “diamagnetic part” of the current, 
— ea(k + eA)I. However, it is proportional to the small pa¬ 
rameter a, and we will show that the main part of the Meissner 
kernel in the superconducting state comes not from this “dia¬ 
magnetic part” but from the inter-band contribution as dis¬ 
cussed in Sect. 3. 

We can carry out a similar calculation of the Meissner ker¬ 
nel as in Sect. 3. Note that the unitary matrix in Eq. (3) does 
not change. After some algebra, we obtain 


C tra ’ (s) ( 0 , 0 ) = 2 e 2 £ 

k,r)=± 


{y 2 k x + r\ak x e{k)) 2 df(E n (k)) 
e 2 (k) dE ri (k) ’ 

(A-3) 


and 


A'™ ter,(s) (0,0) = —2eV£ ( 1 — 

k 

1 


v 2 k 2 

e 2 (k) 


E + (k) +E-(k) 
S £+ lntel (k,A(k)), 

k 


1 - 


- A 2 {k) 


E + (k)E-(k) 


(A-4) 


at T = 0 where E±(k) = ^/^±(k ) 2 + A(k) 2 , and t,±(k) is 
now changed to |±(fc) = ±e(k) + fk 2 — jj.. 

Let us consider the normal state with /.i > M. There are 
several contributions to the Meissner kernel. First is the con¬ 
tribution from the diamagnetic current which appeared in 
Eq. (A-2) artificially. Since this diamagnetic contribution is 
proportional to the electron density, we obtain 


^-dia_^-dia 


K' 


dia 


K dla = 


3n 2 


Mb 


b ia = 


3 n 2 F ' 
(A-5) 

where ko is the artificial Fermi surface of the lower band 
shown in Fig. A l, i.e., ^_(ko) = 0 is satisfied, and kp is the 
Fermi surface of the upper band. Sign ± of the subscript indi¬ 
cates the upper (77 = +) and lower (77 = —) band. The second 
contribution is from the intra-band (Fermi-surface) contribu¬ 
tion from the artificial Fermi surface at ko (Eq. (A-3) with 
A = 0), which is given by 

" 2 ,2(v 2 ko-ak 0 e(ko)) 2 


K ^m m = _±_ K . 


1 


_ _ L’’ 

3 K 2 ' 


h) 


e 2 0o) 


e{k Q ) 


| ako — 


v 2 ko I 
£(*o) 1 


(A- 6 ) 


Note that df/dE- becomes a delta function —S(£-(k)) 

2 

when A = 0, and a — > 0 holds since e(ko) is approxi¬ 

mated as e(ko) ~ vko ~ 2 v 2 /a. Similarly, from the Fermi sur¬ 
face at kp, we obtain 


K 


intra,(N) 


-xx-\- 


m = ( a + 


e(k F ) 


(A-7) 


Finally, the inter-band contribution (0,0) (Eq. (A-4)) 

can be calculated analytically. The integrand in Eq. (A-4) is 
nonzero in the region where ^ + (k)^_(k) < 0 , i.e., in the re¬ 
gion of kp < k(< ko). Thus we obtain 


ip 


-inter,(N) 


(0,0) = £b 

k 


r (M) 



Fig. A-l. (Color online) The dispersion relation of the 4x4 Dirac Hamil¬ 
tonian in normal state with the quadratic term. 
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2n 1 J kF e{k)\ l 3 e(k) 2 ) 

e 2 v 2 k 3 k ° 
kF 


e 1 / v 2 kl v 2 k 3 F \ 
37i 2 \e(k 0 ) e(k F )J 


(A-8) 


We can show that the total of Eqs. (A-5)-(A-8) vanishes, 
which means that the Meissner kernel in the normal state van¬ 
ishes. Next, let us consider the case of superconductivity. We 
assume that the superconducting gap A(k) is finite only for 
k <k c < ko, as discussed in the text. In this case, the artificial 
Fermi surface at ko survives. Thus, the contributions of K d,Ll 
and ^^’(^(OjO) remain. On the other hand, A'^ a '( Ni (0,0) 
vanishes because of the opening of a gap near k = k F . There¬ 
fore, the total Meissner kernel, K, becomes 


K=K^ ev ' {s \o,o)+K dia 


+ K dia + K™ tra ’ {N) (0,0) 


—A, 


■inter, (S) 


( 0 , 0 ) 


e 2 \’ 2 kl ae 2 3 

3^-7(ko) + 3i? kp ' 


(A-9) 


where we have used the explicit results in Eqs. (A-5)-(A-8). 
Equation (A-9) means that we should include the second and 
third terms to obtain the correct Meissner kernel. Since A (k) 
vanishes for k > k c , we can rewrite A'” ter ^ S> (0,0) [Eq. (A-4)] 
as follows. 


4r r ’ (S) (0,0)= £ F inter (k,A(k)) + Y F^(k,A(k)) 

\k\<k c \k\>k c 


= Y F inter (k,A(k)) 
lk\<k c 


e 2 r V 2 k 3 0 

3 n 2 _e(k 0 ) 


v 2 k 3 ' 

W)\ ’ 

(A-10) 


where we have used the similar calculation as in (A-8) for the 
integral of |£| > k c . Then, the total Meissner kernel (A-9) can 
be rewritten as 


K= Y F mtel (k,A(k)) 

\k\<k c 


e 2 v 2 k 3 ae 2 3 

+ 3 n 2 e{k c ) + 3 n 2 F 


(A-11) 


Since a is a small parameter, the last term in (A ll) can be¬ 
come arbitrarily small. Furthermore, it is easily seen that the 
second term is equal to —Y,\k\<k c F mer (k,0), i.e., the normal 
state contribution in the region of 0 < 1*1 < k c . As a result, in 
the small a limit, we obtain 

K = Y E inter (jfc,A(jfc))- Y F intel (k,0), (A-12) 

\k\<k c \k\<k c 


which is equivalent to Eq. (26). 

Although we have used a specific model in this Appendix, 
we think that the subtraction of the normal state contribution 
as in Eq. (26) is reasonable. 
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